Dynamical origin and the pole structure of X(3872) 
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The dynamical mechanism of channel coupling with the decay channels is applied to the case 
of coupled charmonium - DD* states with J PC = 1 ++ . A pole analysis is done and the DD* 
production cross section is calculated in qualitative agreement with experiment. The sharp peak at 
the DqDq threshold and flat background are shown to be due to Breit-Wigner resonance, shifted by 
channel coupling from the original position of 3954 MeV for the 2 s Pi, QQ state. A similar analysis, 
applied to the n = 2, 3 P 2 , 1 Pu 3 Po, allows us to associate the first one with the observed Z(3930) 
J — 2 and explains the destiny of 3 Pq. 



The resonance X (3872) found in Q and confirmed 
and further studied in 0-0] (see @ for review) is still a 
mysterious phenomenon. The measured quantum num- 
bers of X(3872) 01 suggest that this is a 1 ++ state. 
One can list several properties of this resonance which 
are difficult to explain. (1) The width of the peak at 
3872 MeV is zero within experimental energy resolu- 
tion. (2) The peak is exactly at the DqDq threshold 
(3871.2 MeV) and not at a little higher D+D*_ thresh- 
old (3879 MeV); however, isospin conservation predicts 
that both thresholds should enter with the same weight. 

(3) The single-channel theory [tJ predicts a standard 
2 3 Pi level of the QQ system around 3950 MeV; however, 
among the structures observed by Belle in this region, 
X(3940),F(3940),Z(3930), there seems to be no exam- 
ples suggesting the 1 ++ identification @. (4) Among 
the four members of the 2 3,1 Pj multiplct, only one with 
J = 2 can be associated with Z(3930), which was ob- 
served as a regular resonance, A" (3872) looks like a sharp 
cusp, and two others are not seen in this region. It is 
important to explain this very different behavior. 

On the theoretical side there are models based on 
the D D* Q molccularpicture of A(3872) [flUH and the 
tetraquark system 12 1; see [l3[ for a review. However, 
one cannot get a simultaneous explanation of points (1)- 

(4) from these models, and we develop here an alternative 
approach. It is a purpose of this Letter to exploit a re- 



alistic dynamical mechanism, constructed in |14| . which 
can explain all four points. Below we shall Briefly ex- 
plain the mechanism of channel coupling (CC) with the 
decay channels [14|. This method allows us to calculate 



not only position of poles in the CC system, like QQ and 
Qq, Qq, but also scattering amplitudes and production 
cross sections. We demonstrate that, in the 2 3 Pi CC sys- 
tem coupled by the S-wave decays, two poles originating 
from complex conjugate Breit-Wigner resonances of QQ 
system are shifted by CC to the final position, with one 
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pole yielding a narrow cusp at one of thresholds, and an- 
other a shifted flat background. We show that when the 
CC coupling increases, the latter pole yields a very broad 
bump, and at the same time the weak threshold cusp at 
the higher threshold D + D*_ goes over into a sharp peak 
at the lower one DqD^. In this way the same value of 
the coupling constant 7 (well within the accuracy limits 
of the universal constant fitted to different charmonium 
and bottomonium states in [Hj]) produces the visible ef- 
fects, compatible with the properties (l)-(4) mentioned 
above. To produce this effect as in the 2 3 Pi state, the 
original single-channel pole should be above and in the 
attraction region of the threshold. For the poles origi- 
nally below threshold, CC shifts poles down, as it hap- 
pens with the 2 3 P 2 pole. The same approach allows us 
to explain the situation with two other poles, 2 3 Po and 
2 1 Pi, as will be discussed below. 

Resonances in coupled channels can exist for 3 different 
reasons [lj|: (a) due to bound states in the QQ channel, 
which are shifted by CC, (b) due to poles in the (Qq)(Qq) 
channel, shifted by CC; (c) due to strong CC alone (even 
if no interaction exists in decoupled channel). The most 
striking feature of the CC resonance is that it approaches 
the threshold at increasing coupling and typically looks 
like a pronounced cusp at the threshold of small width. 
In the realistic physical problem several of these reasons 
can be present at the same time: e.g., the bare state in 
the single-channel charmonium can be shifted by strong 
CC exactly to the threshold. This situation will be dis- 
cussed below and is characteristic for the single-channel 
pole above the S'-wave threshold, where one meets at 
the same time with the QQ, and the strong CC interac- 
tion, which shifts the bare charmonium pole exactly to 
the threshold position and produces a sharp peak. We 
shall quantitatively describe this situation in the CC for- 
malism of [14j, where the only changeable parameter is 
the channel coupling constant 7 being varied around the 
standard value. 

The basis of the CC theory developed in [l4|, HH can 
be shortly formulated in three relations: a) The effective 
string decay Lagrangian of the 3 Pq type for the decay 
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(Q9)(Qq) 



Csd = / IpqMulpg d 4 



(1) 



with M 
cays [l5 



tested in charmonium and bottomonium de- 
M u sa 0.8 GcV where light quark bispinors 
arc treated in the limit of large m c mass as solutions 
of Dirac equations, and this allows us to go over to 
the reduced (2 x 2) form of the decay matrix element, 
Csd -> jjiv c crpvd 4 x, 7 = {mq+ v^ v+£o) « 1.4 (with 

realistic averages of scalar U = or and vector V — — f 
potentials), b) The decay matrix element of the state 
ni of heavy quarkonium QQ to the states ri2 and 713 of 
heavy-light mesons Qq, Qq, 
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Here all wave functions (^qq, ^Qq, ^Qq) refer to the ra- 
dial parts of the corresponding wave functions, while j/Jfjj 
comprises the decay vertex of L s d and all spin-angular 
parts of mesons involved; the list of yffsj f° r the 6 lowest 
states is given in Table VII of [3. In Eq^c = - " q ^ q , 
where the averaged kinetic energies of heavy and light 
quarks in D meson u) q ~ 0.55 GcV, ujq ~ 1.5 GeV are 
taken from [l7j]. c) The CC interaction "potential" V™^ 
between Qq and mesons due to intermediate states of 
bound QQ system, 



We first apply this method to the case of the 1 ++ state 
of QQ, and we confine ourselves to one state 2 3 Px in the 
QQ system so that Eq. Q reduces to 



E - E n = w nn (E) 



(7) 



and w nn is given in ([5]), where J n n 2 n 3 {p) should be cal- 
culated with the wave functions of D and D* mesons (we 
disregard here the possible difference in wave functions of 
D and D± , and of Dq and D ?j_ , and we also disregard all 
other states n?,, 713 beyond DD*). The wave functions of 
all states involved have been calculated by Badalian et al. 
Q, using the relativistic string Hamiltonian fl8| derived 
in the framework of the field correlator method [19| . Here 
only universal input is used: current quark masses m q , 
string tension a and strong coupling a s . These realistic 
w.f. (^qq, 4>Qq, ipQq) have been fitted by a series of os- 
cillator w.f. and in this way both J n n 2 n 3 {p) and w nn (E) 
were numerically obtained. 

To understand the nature of singularities in the energy 
plane, which produce the cusp at the DqDq threshold, 
we find solutions of (JT)). We separate from w(E) the 
square root singularity, while the rest is a slowly varying 
function, which we approximate by w{E t h) 



w{E) £* w{E th ) 
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where E = E* 



k 2 /2M and M = , a / dM ,?* is the re- 



J t h I 'V 1^ c*x iu ^ m d +m d , 

duced mass. Note that w(E t f l ) < 0, and the two pole 
solutions of Eq.([7]) are 
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and the final (Qq)(Qq) Hamiltonian looks like H = 
Hq + V cc , where Ho may contain the direct Qq and 
Qq interaction, which is 0{\/N c ) and we disregard it in 
what follows. The equation for the pole position is 14 1 



det[E - E - w(E)] = 0, (E) mn = E n S mn (4) 
where E n is the mass of the bare states and w nm (E) is 



w nm (E) 



d^p \ ^ Jnn 2 n 3 (p) Jmn 2 n 3 (p) /r\ 



E-E n2n3 (p) 



with a = M 2 \ J(0)| 2 /tt and b = 2M [w{E th ) + {E n -E th )\. 
Starting with small coupling, one has two Breit-Wigner 
poles. With increasing 7 the square root vanishes at 7 = 
7 and the two poles collide; for 7 > 7, both poles move 
apart along imaginary k axis as shown in FigfT] and at 

sr^rr, the pole k+ passes zero, 



some 7 



,w(E th )\ ' u "° n + 

providing the sharp peak at the threshold. 

The analysis can be extended to the case of two thresh- 
olds eQ and E^ with the resulting equation 



k\ — b + ik\ — 



(K-A) - =0 



with E n2 n 3 (p) = E n2 + E n3 ■ We shall be interested be- 
low in the pole positions, i.e., solutions of (j4|, and the 
most important for comparison with experiment, is the 
production cross section 



Vprod = 2J ^n(E) 



Im w„ n (E) 



\E-E n -Wnn(E)\ 2 ' 



where $„(£') is a weakly changing function of E. In FigsJ5] 
and FigJ3]this factor is omitted. 



2M(£ f ( ( 2) 



E^V). For small A the analysis 



where A = „,„ K ^ th - ^ th 

goes as before, and the only difference in the k\ plane is 
the appearance of the cut connecting points k\ = ±-\/A, 
which denote access to the second sheet of k^- As before, 
the trajectory of the highest pole (k[ + ^) passes through 
(6) the origin, leading to a sharp cusp at E^ . We have 



(2) 
th ' 



found that the pole never passes through the point E 1 

(2) 

implying that the peak at E\ h is never so high, as at 
E^V , compare curves (2) and (4) in FigJ5jD. These curves 
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FIG. 1: Motion of poles fc+ (starting on right-hand side) and 
k- (starting on left-hand side) in the k plane (in units GeV) 
with growing coupling parameter 7. Numbers at field circles 
on curves correspond to numbers on curves in Fig. [2^,. 



correspond to the situations when the pole is closest to 
(2) 

E\ h and when the pole passes the origin, respectively. 

Summarizing this analysis, we have found the pole 
structure behind the phenomenon of V(3872), and we 
may assert that the sharp peak at 3872 is due to the 
pole k[ very close to the DoDq threshold, which origi- 
nally was a genuine Breit- Winger pole, and the flat back- 
ground contains the far virtual pole k[ \ originally the 
complex conjugated Brcit-Wigncr pole generated by the 
same charmonium 2 3 Pi state at ~ 3950 MeV, and shifted 
to the final position k[ ' by CC. We stress that the nec- 
essary condition for the threshold cusp of the type of 
V(3872) is \E„ - E th \ « \w(E th )\, which means that the 
strength of CC should be as large as the distance of the 
original QQ bound state from the threshold. In other 
words, the QQ pole should be "within reach" of CC in- 
teraction. One can also see from (j9|) that for E n < E t h 
the radicand is negative, and for the growing 7 the pole 
k + moves up, farther from threshold, so that for moder- 
ate 7 both poles are far from threshold. 

In the case of two distinct thresholds w(E) contains 
two isotopically equivalent thresholds, which we take into 
account with equal weights. To compare with experi- 
ment, we have used the QQ production cross section ([6]), 
where QQ (in our case cc) are produced in some primary 
reaction, e.g. in e + e~ double charmonium production or 
from B — > KX and then the QQ — > (Qq)(Qq) transition 
takes place. The resulting form of the cross section is 
shown in Figj2]D for five different values of 7. One can 
see from Figj2|b) that for weak coupling [curve (1)] only 
the single-channel charmonium state E n (2 3 Pi) is seen 
with r ~ 35 MeV, and the next two curves display a 
cusp at E th and an almost disappeared QQ resonance, 

while curve ((4]) clearly signals a strong cusp at E$ and 
no other features. At even stronger CC [curve (5)] the 
CC pole goes away from thresholds and the whole picture 
flattens. Thus we see that the experimental situation is 
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(a) One threshold, E th (D Do)=3.872 GeV. 




(b) Two thresholds, E th (D D%; D+D*_ )=3.872; 3.879 GeV. 

FIG. 2: Production [see Eq. © (in units GeV -1 ) for 1++ 
state with different values of channel coupling parameter [(1) 
7 = 0.6, (2) 7 = 1.0, (3) 7 = 1.1, (4) 7 = 1.2, (5) 7 = 1-3]. 
For small values of channel coupling parameter 7 [curve (1)] 
one can see a good Breit- Wigner shape, which corresponds to 
the shifted 2 3 Pi state, while for larger 7 [curves (3) and (4)] 
there is a broadening of higher resonance together with steep 
rise near the threshold E th (D Do) = 3.872 GeV. 



well reproduced by curve (4) . The positions of both poles 
changing with 7 are marked in Fig[TJ One can easily see 
how the pole k + produces the sharp cusp in position (3) 
for one threshold treatment, corresponding to curve (3) 
in the production cross section in FigJ^a). 

Having found the mechanism, producing the peak at 
the DqDq threshold, one may wonder what happens with 
other states of the n = 2 3 Pj family, J = 0, 2. To this 
end one should first estimate the position of bare poles 
EL (see Table IJ). We use the results of Badalian et al. 
[Tj with the slightly modified spin-orbit interaction. 

We take now the 2 3 P2 bare state, which is mostly con- 
nected with the D*D* channel, while the DD state is 
in the D wave and can be neglected. One can see, that 
E n (2 3 P 2 ) < E th (D*D*), and w(E) in © is real and neg- 
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TABLE I: Hadronic shift (MeV) of charmonium 2 J ' i P J bare 
states E n for different channels. The bare positions were 
taken from Badalian et al. [?|, 5 is the total shift and n 
is the closed channel (D*D*) suppression coefficient. 



State 






Shifts (7 = 1.1) 
DD DD* D*D* 8 


E Exp. 


2 3 P 2 


2++ 


3969 


- K = 

- k = 0.25 -14 

- k = 0.5 -27 




-14 
-27 


3.969 

3.955 Z(3930) 
3.942 


2 3 Pi 


1++ 


3954 


to threshold 


3.872 X(3872) 


2 3 P 


0++ 


3918 


-2 - k = 
-2 - k = 0.25 -33 
-2 - k = 0.5 -66 


-2 
-35 
-68 


3.916 
3.883 
3.850 


2 1 P 1 


1+- 


3959 


- -25 k = 

- -29 n = 0.25 -7 

- -32 k = 0.5 -14 


-25 
-36 
-46 


3.934 
3.927 
3.918 



ative. Hence, with increasing coupling the pole is shifted 
down, away from the threshold. This is shown in Ta- 
ble |TJ In terms of our previous analysis, using Eq. J9|) 
with b < 0, one can see that the pole is on the imaginary 
axis. Moreover, one can estimate that y/\b\ 3> § and 
the near-the-threshold approximation ([5]) is not applica- 
ble; one can better use the original Eq. ([7|) , which yields 
a shift of |?£)(.E n )| ~ 55 MeV. At this point one should 
take into account the necessity of renormalizing the con- 
tributions of higher closed thresholds (which otherwise 
produce unacceptable shifts, see [3,[20, Therefore 
we introduce the coefficient k, which multiplies in (j3|) the 
contribution of the closed channel D*D*. We estimate 
K in an approximate range < n < 0.5 because it gives 
the most sensible results for the mass shift of J/ip. How- 
ever k may depend on quantum state and bare energy. 
The resulting position is near the experimentally found 
@ Z(3930) peak, as seen in Table Q] 

A similar situation occurs for the 2 3 Pq state, where the 
lower threshold is far below, |P„(2 3 P )-£: t/l (Z)P>)| - 200 
MeV, while the higher threshold is more distant, than in 
the case of the 2 3 P2 state. Our calculation for the shift 
of the 2 3 P state yields a large value AE ~ 68 MeV for 
K = 0.5, see Tabled so the final position is 3.850 MeV, 
which possibly corresponds to the position of the wide 



peak in e 



J/i(jDD in [211 ]. This enhancement of 



the shift is due to a much larger overlap matrix element 
in the 3 Po state. A similar situation was discussed in 

An interesting situation occurs in the case of the 1 + 
state with the 2 1 P\ pole at bare position 3959 MeV. Here 
the coupling to the DD* channel is much weaker, than 
in the case of the 1 ++ state, and the main coupling is 
with the closed D*D* channel, which defines the destiny 
of the bare pole. When we renormalize the contribution 
of the D*D* channel with the coefficient discussed 
above, the final shift for k = 0.5 is around 45 MeV. In 
Fig. [3] we demonstrate how the production cross section 
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FIG. 3: Production (see Eq. [in units GeV -1 ] for the l + ~ 
state in DD* channel for channel coupling parameter 7 = 1.1 
and with different values of closed D* D* channel suppression 
coefficient, n — 0; 0.25; 0.5. 



changes with k, and one can see, that the resulting width 
at k = 0.25 ~ 0.5 is around T ~ 50 MeV. 

In conclusion, we have calculated the amplitudes of 
CC processes connecting QQ and (Qq)(Qq) systems via 
the decay matrix element, Eq.([3]), involving realistic wave 
functions of all hadrons involved, and the CC constant 
7(M W ), fitted earlier to bottomonium and charmonium 
transitions. For the concrete case of the 1 ++ state of 
charmonium we have found pole structure and produc- 
tion cross section. At small CC two poles correspond to 
the complex conjugated poles of one Breit- Winger reso- 
nance of the 2 3 Pi state of QQ with the width r ~ 35 
MeV. For increasing CC one of the poles approaches the 
thresholds and another moves away, as a result this bare 
resonance flattens, while a sharp cusp appears first at the 
D+D*_ and then at the D D£ threshold at - 3872 MeV. 
This latter situation with the sharp narrow cusp at the 
D n D*. threshold and absence of any other structures (ex- 
cept for a tiny cusp at higher threshold) including the 
region around 3940 MeV corresponds to the observed 
production yield @. We conclude that our dynamical 
mechanism explains properties (l)-(4), in particular, why 
the resonance V(3872) is at the lower, but not the higher 
threshold, why it is so narrow, and why the original 2 3 Pi 
state of charmonium is not seen in experiment. 

An alternative and close in spirit approach was dcvcl- 
H, [H|. 



oped recently in 
the conclusion in 
cc 2 3 Pi state origin' 

m, 



Our analysis partly supports 
24(, that "A" (3872) may be of ordinary 
Our results differ from those of 
where two 1 ++ states were found, one associated 
with X (3872), and another with A" (3940). In a recent 
review ^ the CC analysis of the V(3872) and V(3940) 
was reported with the conclusion, that both states can- 
not be reproduced in the exploited model simultaneously. 
This result is in common with ours, since in our case the 
broad enhancement due to the second pole is near the 
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DqDq threshold and cannot be associated with Y(3940). 
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